High-order harmonic generation by enhanced plasmonic near-fields in metal 

nanoparticules 
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We present theoretical investigations of high-order harmonic generation (HHG) resulting from 
the interaction of noble gases with localized surface plasmons. These plasmonic fields are produced 
when a metal nanoparticle is subject to a few-cycle laser pulse. The enhanced field, which largely 
depends on the geometrical shape of the metallic structure, has a strong spatial dependency. We 
demonstrate that the strong non-homogeneity of this laser field plays an important role in the HHG 
process and leads to a significant increase of the harmonic cut-off energy. In order to understand 
and characterize this new feature, we include the functional form of the laser electric field obtained 
from recent attosecond streaking experiments [F. Siifimann and M. F. Kling, Proc. of SPIE, Vol. 
8096, 80961C (201 1)] in the time dependent Schrodinger equation (TDSE). By performing classical 
simulations of the HHG process we show consistency between them and the quantum mechanical 
predictions. These allow us to understand the origin of the extended harmonic spectra as a selection 
of particular trajectory sets. The use of metal nanoparticles shall pave a completely new way of 
generating coherent XUV light with a laser field which characteristics can be synthesized locally. 

PACS numbers: 42.65.Ky,78.67.Bf, 32.80.Rm 



When matter, i.e. atoms or molecules, is exposed to 
short and intense laser radiation, non-linear phenom- 
ena are triggered as a consequence of this interaction. 
Amongst these phenomena, high-order harmonics gen- 
eration (HHG) process [l], Q has attracted considerable 
interests, since it is one of the most reliable pathways 
to generate coherent light from the ultraviolet (UV) to 
extreme ultraviolet (XUV) spectral range. As a result, 
HHG has proven to be a robust source for the generation 
of a PHz attosecond pulses train 0], that can be tem- 
porally confined to a single XUV attosecond pulse, now 
with kHz repetition rates [J- Thanks to its remarkable 
properties, HHG can be used as well to extract temporal 
and spatial information with both attosecond and sub- 
Angstrom resolution on the generating system Hence, 
HHG represents a considerable tool to enable scrutinizing 
the atomic world with its natural temporal and spatial 
scales 

The intuitive physical mechanism behind HHG, for a 
single atom or molecule (referred to as 'single emitter'), 
has been well established in the so-called three steps or 
simple man's model p^4I3 |: in the first step, an elec- 
tronic wave packet is released the continuum by tunnel 
ionization through the potential barrier as a consequence 
of the non-perturbative interaction of the single emitter 
with the laser field. In the second step, the emitted elec- 
tronic wave packet propagates away from its ionic core in 
the continuum to be finally driven back when the laser 
electric field changes its sign. In the final step, upon its 



return, the electronic wave packet may recombine with 
the core and the system relaxes the excess kinetic energy 
acquired by radiating a high-harmonic photon. 

In order to experimentally control the high harmonic 
features two main types of approaches have been at- 
tempted. The first one is based on the control of the 
HHG process via the manipulation of the laser field char- 
acteristics in time and/or space [l5[. The second one is 
based on the control of macroscopic properties of the tar- 
get samples (i.e. phase matching) leading to very inge- 
nious target geometries fl6Hl"8l |. However, all these ap- 
proaches rely on detecting the far-field properties of the 
harmonic yield which is a consequence of collective single 
emitter's ones. Therefore, it is legitimate to study how 
a synthesized single emitter could lead to new parame- 
ters for HHG control. Metal nanoparticles are a matter 
of choice since the spatial geometry of the nanoparticle 
and the material used can be chosen to confer a spatial 
transverse non-homogeneity to the laser field. One of 
the first demonstrations of such an effect was obtained 
for surface plasmonic resonances that can locally am- 
plify the laser field [lj| . The local electric fields resulting 
from such resonances can reach an enhancement greater 
than 20 dB (2^, HH . Consequently, when a femtosecond 
low intensity laser pulse couples to the plasmonic mode 
of the metal nanoparticle, it initiates a collective oscil- 
lation among free charges (essentially electrons) within 
the metal. A location of highly amplified electric field 
is thus created while these free charges redistribute this 
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field around the metal nanostructure. The enhanced field 
is well above the threshold for generating high harmonics. 
So, by injecting noble gases surrounding the nanoparticlc, 
HHG can be produced. Particularly, whilst using gold 
bow-tie shaped nanostructures, it has been demonstrated 
that the initial laser field (800 nm laser with intensity 
10 11 W/cm 2 ) can be enhanced sufficiently to generate 
XUV photons, i.e. from the 7th (114 nm) to the 21st (38 
nm) harmonic order. Furthermore, the high harmonics 
radiation generated from enhanced laser field, localized 
at each nanostructure, acts as a point-like source, en- 
abling collimation or focusing of this coherent radiation 
by means of constructive interference. This localization 
of the enhancement confers to the transverse laser field 
profile its non-homogeneity in the region where the elec- 
tron dynamics is taking place. Besides, the interaction 
length along the propagation direction of the laser field is 
restricted to few nanometers so that no phase conditions 
needs to be considered in order to observe and calculate 
the emitted harmonic yield. As a consequence, spatially 
arranged nanostructures open a wide range of possibili- 
ties to enhance or shape the spectral and spatial prop- 
erties of the fundamental laser field and the harmonic 
field These two features imply strong modifications 
in the harmonic spectra which will raise the interest of 
the strong field community to utilize nanoparticles sur- 
rounded by gas atoms or molecules as a new type of tar- 
get [2214251 ]. However, the initial thrill about the utiliza- 
tion of plasmonic fields for HHG in the XUV range, was 
put in debate by recent findings [26l - |28j . Fortunately, al- 
ternative ways to amplify coherent light using plasmons 
were explored. To cite only a couple of prominent exam- 
ples, the production of high energy photoelectrons using 
plasmonic near enhanced fields from dielectric nanopar- 
ticles I2II, m etal nanoparticles [30M32| | and metal nan- 
otips (33p[37 1 appears to be perfectly plausible. Besides 
the question of damage threshold of such nanotargcts has 
been highlighted and indicates that new routes using ini- 
tially low intensity laser field need to be considered [38| . 

In this contribution we investigate how HHG yield can 
be controlled using enhanced near-fields. These fields are 
obtained only when a metal spherical nanoparticlc is il- 
luminated by a few-cycle laser field . This particularly 
ultra-short interaction regime confers to these near-fields 
a strong spatial dependence that can be accessible ex- 
perimentally using attosecond streaking techniques [39| . 
The plasmonic enhanced near- field in the vicinity of these 
nanotargcts could act as a femtosecond coherent source 
and therefore drive strong laser-matter processes in the 
surrounding atomic or molecular gas , e.g. high-order 
harmonic generation. In addition, the enhancement of 
the fundamental laser field reaching one or more orders 
of magnitude, depending on the size and material em- 
ployed, would allows us to consider low input laser in- 
tensity and thus to work well below the damage thresh- 
old of the nanotargct. Using both quantum mechanical 
and classical approaches, we predict that the signature 
of such synthesized fields is pronounced in the harmon- 



ics cutoffs which now extend far beyond the conventional 
semiclassical limits. 

Most of the numerical and semiclassical approaches to 
study high-order harmonic generation (HHG) are largely 
based on the assumption that the laser electric field 
(E(r, t)) and its associated vector potential (A(r,i)) are 
spatially homogeneous in the region where the electron 
dynamics takes place, i.e. E(r, t) = E(£) and A(r,t) = 
A(t) [10, Ej]. Nonetheless, near-fields generated in the 
vicinity of metal nanoparticles are not spatially homo- 
geneous and we address the question on how such non- 
homogeneity could be revealed in the harmonic spectra. 
In the case of homogeneous fields, the HHG process has 
been theoretically tackled using different approaches (for 
a summary see e.g. [42|, |43[ and references therein). In 
this article, we compute the HHG spectra by including 
the spatial dependence of the field in the dimcnsionally 
reduced Time Dependent Schrodinger Equation (TDSE) 
by considering the actual functional form of the laser 
electric field spatial dependence, obtained from attosec- 
ond streaking experiments. The TDSE in one spatial 
dimension and for a model atom can be written as 1401: 
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In here, V a {x) is the atomic potential and Vi(x,t) repre- 
sents the potential due to the laser electric field. For the 
V a (x) , we use a soft-core potential to avoid singularities. 



V a (x) = - 



1 



Vx 2 + 1 



(2) 



This potential was first proposed in [44( and has been 
widely used in studies of laser-matter processes in atoms. 
Certainly the model potential is not suitable to predict 
the actual atomic structural information present in the 
harmonic spectrum, but it is completely applicable to 
characterize the HHG cutoff, once the ionization poten- 
tial of a given atom is set. A particular ionization poten- 
tial can be defined by varying the parameter a in Eq. © . 
In this work we use a = 1.62 to model Xenon atoms with 
Ip = 12.1299 eV (0.446 a.u.). In addition, we assume that 
the xenon atom is in its ground state before the laser field 
(t = —00) is turned on. Equation ([1]) is solved numeri- 
cally by using the Crank-Nicolson method [40| • To avoid 
spurious reflections from the spatial boundaries, at each 
time step, the electron wave function is multiplied by a 
mask function [45[ . The potential due to the linearly po- 
larized laser electric field in the a;-axis, E(x,t) in ([T]), is 
given by: 



Vi(x,t) = -E(x,t)x. 



(3) 



We employ the function given by [30( to define E(x,t), 
i.e. 



E(x, t) = E f(t) exp(— x/x) sm(ut + 



(4) 
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In Eq. (j4]), Eq, uj, f{t) and (j) are the peak ampli- 
tude, the laser field frequency , the field envelope and 
the carrier envelope phase (CEP), respectively. The spa- 
tial dependence of the plasmonic near-field is given by 
the term exp(— x/x) and it depends on both the size and 
the material of the spherical nanoparticle used. Eq. (|4]) 
is valid for x > i?o, where Rq is the radius of the metal 
nanoparticle, i.e. for values of x outside of the metal 
nanoparticle. Additionally, it is important to note that 
the electron motion stands in the region x > Rq with 
(x + R ) ^> 0. In this work we consider the laser field 
having a sin 2 envelope to avoid DC components: 



tot 



(5) 



where n p is the total number of optical cycles so that the 
total pulse duration is r = 2im p /uj. 

The harmonic yield of the atom is obtained by Fourier 
transforming the acceleration a(t) of the electronic wave 
packet HI]: 



D(w) 



1 1 



dtc~ iult a(t) 



(6) 




with a(t) obtained by using the commutator relation: 
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FIG. 1: High-order harmonic generation (HHG) spectra for 
Xe (ionization potential I p = —0.446 a.u.), laser wavelength 
A = 720 nm and intensity I = 2 x 10 13 W-cm~ 2 . We use a sin 2 
shaped pulse with n = 5 (about 13 fs of total time and « 6.5 
fs FWHM). Panel (a) shows the homogeneous case, panel (b) 
X = 50 and panel (c) \ = 40. The arrow in panel (a) indicates 
the cutoff predicted by the semiclassical model [13 ] . 
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In here, %{t) and &(x, t) are the Hamiltonian and the 
electron wave function defined in Eq. ([T]), respectively. 
The function D(u>) is known as the dipole spectrum, 
which gives rise to the spectral profile measured in HHG 
experiments. 

Figure 1 depicts the harmonic spectra for Xenon gen- 
erated by a laser pulse with I = 2 x 10 13 W/cm 2 , wave- 
length A = 720 nm and a total pulse duration of 13 fs, 
i.e rip = 5 (this pulse corresponds to « 6.5 fs FWHM). 
We have chosen these parameters to perform calculations 
consistently with experimental conditions (30j . For the 
case of a homogeneous field, with this low fundamen- 
tal laser intensity, no harmonic beyond the 9 th order 
are observed (see the arrow in panel a). The spatial 
decay constant \ quantifies the non-homogeneity due 
to the nanoparticle. It varies together with the size of 
the particle and the metal employed. Changing \ is 
therefore equivalent to choosing the type of nanopar- 
ticle used, which leads to overcome the semiclassically 
predicted cutoff limit and reach higher harmonic orders. 
For instance with \ = 40, harmonics in the mid 20s are 
obtainable (panel c) and with \ — 50 we observe har- 
monic orders well above the 9 th (a clear cutoff at n « 15 
is visible) (panel b). In addition to this cut-off exten- 
sion, we observe a change in the harmonic periodicity. 
This is related to the breaking of symmetry imposed by 
the induced non-homogeneity of laser electric field (see. 
Eq. O). 

In order to understand the relationship between the 
laser field induced non-homogeneity and the harmonic 



cutoff extension, we consider the semi-classical three- step 
model [HI, IHi • As it was already pointed out in [U, |47j . 
this new feature may appear due to the combination of 
two main factors, namely the non-homogeneous charac- 
ter of the laser electric field and the confinement of the 
electronic motion. For an homogeneous field, it is well 
established that the position of the high-order harmonic 
cutoff holds as: 



n c = (3.17f/ p + I p )/u> 



(8) 



where n c is the harmonic order at the cutoff, u> the laser 
frequency, U p the pondcromotive energy (defined by U p = 
I/Aoj 2 , with I being the laser intensity in a.u.) and I p 
the ionization potential of the atom or molecule. In our 
case we solve numerically the Newton equation for an 
electron moving in a linearly polarized (in the a;-axis) 
electric field with the same parameters used in the 1D- 
TDSE calculations, i.e. with a time and space dependent 
electric field E(x,t) of the form Eq. Then we find 
the numerical solution of: 



x(t) = S7 x Vi aser (x,t) 

= -E{x,t) - [V x E{x,t)]x{t) 

X 



(9) 



where E(x,t) is defined by Eq. (gj. Note that in ©, 
we have presented explicitly the time dependence of the 
electron trajectory, i.e. x = x(t). We consider initial con- 
ditions similar to the simple man's model: the electron 
starts at position zero at t — ti (the ionization time) with 
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FIG. 2: (Color online) Total energy of the electron (expressed 
in harmonic order) driven by the laser field calculated from 
Newton's second law and plotted as a function of the ioniza- 
tion time (green circles) or the recollision time (red circles), 
(a) Homogeneous case, (b) x = 50 and (c) x — 40. 



zero velocity, i.e. x{tj) = and i(to) = 0. When the elec- 
tric field reverses its direction, the electron returns to its 
initial position (i.e the electron recollides or recombines 
with the parent ion) at a later time t = t r (the recollision 
time), i.e. x(t r ) = 0. The electron kinetic energy at the 
return time t r is calculated from: 



E k (t r ) 



±(t r y 



(10) 



and finding the value of t r (as a function of U) that max- 
imizes this energy, fulfills Eq. ([8]). 

Fixing the value of the ionization time U it is possi- 
ble to compute the classical trajectories and to numeri- 
cally calculate the recollision times t r , i.e. the t r when 
x(t r ) = 0. For a given ionization time U the electron 
trajectory is completely determined, because the second 
order differential equation ([9]) and its initial conditions 
have a unique solution. In Fig. 2 panels (a)-(c), we show 
the dependence of the harmonic order upon the ioniza- 
tion time (ti) and recollision time (t r ), calculated from 
n = {E k {t it r) + I p )/oj with I = 2x 10 13 W cm~ 2 , A = 720 
nm and a sin-squared shaped pulse of 5 cycles (total time 
duration 13 fs). 

Panels (a), (b) and (c) depicts the cases of (\ —> oo) 
(homogeneous field), \ = 50 and \ = 40 (two cases of 
non-homogeneous field), respectively. From panel (a), 



we observe that the maximum kinetic energy of the re- 
turning electron agrees with Eq. (J5J (no harmonic order 
beyond n c ~ 9lj is reached). On the other hand, panels 
(b) and (c) show how the non-homogeneity of the laser 
field modifies considerably the electron trajectories to- 
wards an extension of the harmonic cut-off energy. This 
is clearly present at n c ~ 18cj (28 eV) and n c ~ 27a; 
(42 eV) for x = 50 and \ = 40, respectively These last 
two cutoff values are indeed consistent with the quantum 
mechanical calculations presented in Fig. 1 

In order to show how our model behaves for higher 
laser intensities we compute the harmonic spectra by in- 
creasing the intensity to I = 5 x 10 13 W/cm 2 , keeping all 
other parameters the same (we note that the saturation 
intensity of xenon is w 8 x 10 13 W/cm 2 and consequently 
our proposed values are well below this value). The re- 
sults are shown in Fig. 3. In panels (b) and (c), the non- 
homogeneity of the laser field manifests itself as a clear 
harmonic cutoff extension, reaching values of n c ~ 60w 
(93 eV in energy) for \ = 40. These values show a 
highly nonlinear dependence of the harmonic cutoff with 
the spatial decay constant x- This non-linear behavior 
can be exploited to generate high order harmonic in the 
XUV regime using modest laser intensities. Finally, in 
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FIG. 3: Idem Fig. 1 but now the laser intensity is I = 5 x 10 
W/cm 2 . 

Fig. 4, we plot the recollision time t r of the electron as 
a function of the ionization time U as for the case pre- 
sented in Fig. 1. For the ionization times confined be- 
tween 0.5 < t t < 1.5 optical cycles, the long trajectories 
are those with recollision times t r > 2.25 optical cycles 
and they are visible only for the homogeneous case (red 
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squares (■)) and for \ = 50 (green filled circles (•)). On 
the other hand, short trajectories are characterized by 
t r < 2.25 optical cycles and these are present for both 
the homogeneous and non-homogeneous cases. Our re- 
sults are consistent with those presented in [23|, [24|, [4?J , 
although in our work we use a different functional form of 
the electric field resulting from experimental results. In 
Fig. 4, we observe how the long trajectories are strongly 
modified by the non-homogeneity of the laser field. In- 
deed the homogeneous long trajectories (red squares (H)) 
with ionization times tj around the 1.25 and 1.75 optical 
cycles converge with the short one into a unique trajec- 
tory set (A). Moreover, the branch with tj ~ 1.75, has 
now ionization times smaller than in the homogeneous 
case, hence, the propagation time of the electron in the 
continuum increases so that it can gain a higher amount 
of kinetic energy [2^, [24|, 113: as confirmed in the classical 
calculations presented in Fig. 2. 

The non-homogeneity of the laser field annihilates the 
long trajectories and only short trajectories are now re- 
sponsible for the harmonic spectrum. This behavior is 
also observed using a non- homogeneous field with spa- 
tial linear dependence (for details see Hi, II E|). 
Considering that the electric field strength at the ioniza- 
tion time for short trajectories is higher than for the long 
ones, and taking into account that the ionization rate is 
a highly nonlinear function of this electric field [4!| [5(| , 
long trajectories are much less efficient than the short 
ones. Accordingly, this fact explains why the harmonic 
spectrum, obtained with a non-homogeneous laser field 
that modified the trajectories shows an extended cutoff. 



In conclusion, we present how the high-order harmonic 
generation from xenon atoms is modified by using a 
plasmonic near enhanced field generated when a metal 
nanoparticle is illuminated by a short and intense laser 
pulse. The functional form of the resulting laser electric 
field is extracted from attosecond streaking experiments 
and incorporated in our quantum mechanical approach. 
We employ the numerical solution of the time dependent 
Schrodinger equation (TDSE) in reduced dimensions in 
order to calculate the harmonic spectra. We observe 
an extension in the harmonic cutoff position that could 
lead to the production of XUV coherent laser sources 
and opening the avenue to the generation of attosecond 
pulses from spatially inhomogeneous laser fields. This 
new feature is a consequence of the induced laser field 
non-homogeneity only, which modifies substantially the 
electron trajectories. Furthermore, our quantum me- 
chanical numerical results agree with the classical sim- 
ulations. A more pronounced increment in the harmonic 
cutoff, in addition with an appreciable growth in the con- 
version efficiency, could be reached varying, for instance, 
the radius of the spherical metal nanoparticles and by 
choosing the adequate metal materials. These new de- 
grees of freedom could pave the way to enhance the har- 
monic spectra reaching the XUV regime with modest in- 
put laser intensities. 




1 2 3 

Ionization time (optical cycles) 

FIG. 4: (Color online) Dependence of the semi-classical tra- 
jectories on the ionization (U) and recollision (t r ) times for 
different values of x- Red squares are homogeneous case, i.e. 
X — ► oo; green circles are \ = 50 and blue triangles \ = 40. 
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